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Abstract. - The energy dependence of charged-hadron production in proton-proton collisions at 
RHIC and LHC energies is investigated in a nonequilibrium-statistical relativistic diffusion model 
(RDM) with three sources for particle production. Calculated charged-hadron pseudorapidity 
distributions for pp at RHIC energies of y/s — 0.2 and 0.41 TeV, and at LHC energies of 0.9, 
2.36 and 7 TeV are optimized with respect to the available data. Predictions for 14 TeV are 
made. The central source arising from gluon-gluon collisions becomes the major origin of particle 
production at LHC energies. The midrapidity dip is essentially determined by the interplay of the 
three sources. 
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Introduction. — The investigation of particle pro- 
duction in proton-proton collisions at LHC energies is ex- 
pected to yield new insights into the underlying partonic 
processes. Data from the experimental collaborations are 
now available starting at the injection energy of s/s = 0.9 
TeV, via 2.36 TeV, to the current maximum energy of 7 
TeV 

A particularly interesting observable is the charged- 
hadron multiplicity density per unit of pseudorapidity. At 
midrapidity, it was found to be about 15% higher than 
predicted by the available Monte Carlo models that had 
been calibrated at Tevatron energies [3j. The distribution 
functions for non-singlc-diffractive events have also been 
measured away from midrapidity, with |?y| < 2.5 achieved 
so far jE]. Their shapes are sensitive to the partonic pro- 
cesses that are responsible for charged-hadron production. 

In this Letter I propose to analyze the pseudorapid- 
ity distribution functions of produced charged hadrons 
in pp— collisions at RHIC and LHC energies as measured 
by the PHOBOS |6 and CMS [IJ[5] collaborations in 
a schematic nonequilibrium-statistical model with three 
sources. 

Similar ALICE data at LHC energies are also available 
0, as well as older UA5 data [8] at 0.9 TeV. Correspond- 
ing ATLAS results [2f3| can not be compared directly with 
the CMS and ALICE data because particles and events are 
selected in different regions of phase space. 



The relativistic diffusion model (RDM) has proven to 
be useful in describing and predicting pseudorapidity dis- 
tributions of produced charged particles in heavy-ion coll- 
sions at SPS, RHIC and LHC energies |9J. Related models 
had also been used in low-energy (non-relativistic) heavy- 
ion physics |10| . 

In heavy-ion collisions, the number of produced charged 
hadrons is much bigger than in pp - of the order of 
20,000 charged hadrons in a central PbPb collision at 
y / ijv r iv= 2. 76 TeV - and consequently, the application of 



nonequilibrium-statistical methods such as 9pT is clearly 
justified. 

Special QCD-effects such as the coherence of soft glu- 
ons 12 that had been predicted by perturbative QCD at 
to produce visible effects in charged-hadron 
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distributions generated by e + e~ collisions are less impor- 
tant in the heavy-ion environment since these are averaged 
out through the random properties of the many-particle 
system. 

Proton-proton collisions at the current maximum LHC 
energy of 7 TeV produce about 70 charged hadrons inte- 
grated over the full rapidity space, including the unmea- 
sured region. Soft-gluon coherence as well as other coher- 
ent QCD-effects may still be visible in the data, although 
less pronounced than in electron-positron collisions. 

The number of produced particles in pp is prob- 
ably already large enough to test the usefulness of 
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noncquilibrium-statistical concepts, although it may be 
difficult to observe many-particle effects such as the shift 
of the fragmentation-peak positions towards midrapidity 
with decreasing cm. energy that is clearly seen in the 
heavy-ion data, and can be described analytically in the 
relativistic diffusion model [9 '. 
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Within the RDM, 1 investigate in this Letter the energy 
dependence of the three sources for particle production 
in proton-proton collisions at RHIC and LHC energies. 
The energy range considered here covers RHIC energies of 
t^snn — 0.2 and 0.41 TeV, the presently accessible LHC 
energies of 0.9, 2.36 and 7 TeV, and the maximum LHC 
energy of 14 TeV. 

The sources correspond to gluon-gluon induced pro- 
duction of charged hadrons centered at midrapidity, and 
quark-gluon processes centered at large rapidities, typi- 
cally < y\ : 2 >— T2.8 at 7 TeV. Their relative sizes deter- 
mine the midrapidity dip in the charged-hadron pseudora- 
pidity distributions when added incoherently. The model 
is complementary to QCD-based approaches that rely on 
the corresponding partonic structure functions. 

Prominent and detailed models for multiple hadron pro- 
duction are available in the literature. In particular, the 
dual parton model (PPM ) (16 17 and the equivalent 
quark-gluon string model |18f|20| are based on the cre- 
ation and breaking of quark-gluon strings. 

There the total inclusive hadron production cross sec- 
tion in pp collisions at energies in and below the RHIC en- 
ergy range arises from contributions of two quark-diquark 
chains that overlap in rapidity space. These yield the total 
inclusive cross section when added up incoherently, with 
a minimum at midrapidity. At higher energies, also mul- 
tichain contributions become significant that are likely to 
contribute in the midrapidity region. There are no trans- 
port effects considered in the model. 

Although there is presently no direct connection to the 
DPM, the 3-sources RDM provides an analytical frame- 
work to investigate the interplay of central and fragmen- 
tation sources, transport effects, and their dependence on 
incident energy. 

The model is considered in Sec. 2, the calculation of 
pseudorapidity distributions of produced charged hadrons 
in Sec. 3, and conclusions are drawn in Sec. 4. 

Linear Relativistic Diffusion Model. The Rel- 
ativistic Diffusion Model (RDM) has been developed to 
deal with ensembles of many particles and their distribu- 
tion functions in transverse momentum and rapidity space 
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In particular, it is well-suited to predict and de- 
scribe charged-hadron rapidity distributions in relativis- 
tic heavy-ion collisions from AGS, SPS and RHIC ener- 
gies, to LHC energies [9j. It is tested here for proton- 
proton collisions at LHC energies, where the number of 
produced charged hadrons appears to be sufficiently large 
for nonequilibrium-statistical concepts to apply. 

The rapidity distribution of produced particles emerges 
from an incoherent superposition of the beam-like frag- 



mentation components at larger rapidities arising mostly 
from valence quark-gluon interactions, and a component 
centered at midrapidity that is essentially due to gluon- 
gluon collisions. All three distributions are broadened in 
rapidity space as a consequence of diffusion-like processes. 

The time evolution of the distribution functions is gov- 
erned by a Fokker-Planck equation (FPE) in rapidity 
space [ll] (and references therein). In the linear model, 
it is formulated as an Uhlenbeck-Ornstein [2l] process, 
applied to the relativistic invariant rapidity for the three 
components Rk(y, t) (fc=l,2,3) of the distribution function 
in rapidity space. For a symmetric system such as pp this 
becomes 



d_ 

Ft 



Rk(y,t) 



d_ 
dy 

dy 2 



J(y)-Rk(y,t) 



D k y -R k {y,t) 



(1) 



with the rapidity y = 0.5 • ln((E + p)/{E - p)). The 
beam rapidity can also be written as ybeam = Tymax — 
=F ln(-y/s/ra p ). The rapidity diffusion coefficient D y that 
contains the microscopic physics accounts for the broad- 
ening of the rapidity distributions. The drift function is 



J(y) = -y/i-y 



(2) 



with the rapidity relaxation time t v . It determines the 
shift of the mean rapidities towards the central value with 
increasing time, 

Since the equation is linear, a superposition of the dis- 
tribution functions [22j 23 using the initial conditions 
Ri,2(y,t = 0) = 5(y ± y ma x) with the absolute value of 
the beam rapidities y ma x, and Rz(y,t = 0) = S(y) yields 
the exact solution for a symmetric system. 

In the solution, the mean values are obtained analyti- 
cally from the moments equations as 



< 2/1,2(0 >= TVrnax exp (~t / T y ) 



(3) 



for the sources (1) and (2) with the absolute value of the 
beam rapidity y m ax- The gluon-gluon source remains cen- 
tered at for pp collisions, or other symmetric systems. 
Both mean values < j/1.2 > would attain y=0 for t— > 00, 
whereas for short times they remain between beam and 
equilibrium values. The variances are 



<2, e «(t) = Dl^r y [l - exp(-2t/r„)] 



(4) 



and the corresponding FWHM-values are obtained from 
T = \/81n2 • a since the partial distribution functions 
are Gaussians in rapidity space (but not in pseudorapidity 
space). 

The midrapidity source that arises from gluon-gluon in- 
teractions with mean value zero comes close to thermal 
equilibrium with respect to the variable rapidity during 
the interaction time Tint] the width approaches equilib- 
rium twice as fast as the mean value. I use the notion 
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R gg (y,t) for the associated partial distribution function 
in y-space, with charged particles, cf. table [l] 

Full equilibrium as determined by the temperature 
would be reached for Ti nt /T y 3> 1. The centers of the 
fragmentation sources would then move to midrapidity ac- 
cording to the solution of the FPE, the incoherent sum of 
the three sources would reach a thermal distribution in 
pseudorapidity space. For finite times, however, the frag- 
mentation sources do not reach < 7/1.2 >= during the 
interaction time and hence, remain far from thermal dis- 
tributions in rapidity space, and do not equilibrate with 
the central source. 

Pseudorapidity distributions. If particle identi- 
fication is not available, one has to convert the results to 
pseudorapidity, r\ = — ln[tan(0/2)] with the scattering an- 
gle 8. The conversion from y— to r\— space of the rapidity 
density 



dN 
drj 



dN dy 
dy drj 



p dN 
E^y~ 



J(r), (m)/{p T )) 



dN 
dy 



(5) 



is performed here through the approximated Jacobian 



J(V, M/(pt)) = cosh(??)- 

[l + {{m)/{p T )f 



sinh 2 (7y)] 



-1/2 



(6) 



The average mass < m > of produced charged hadrons 
in the central region is approximated by the pion mass 
since pions represent by far the largest fraction of pro- 
duced charged hadrons, in particular in the midrapidity 
source where the transformation has the biggest effect. 
The mean mass is larger in the fragmentation region, with 
< m >~ mp/n c ' h + ■ (n^ — l)/« c ^ — 0.27 GeV at 
Js = 7 TeV where n\? = dN/dr>(< m 2 >) ~ 6. 



l lh = dN/dr,(< m,a > 
Due to the Jacobian, the partial distribution functions 
differ from Gaussians, but as a consequence of the rela- 
tively high mean transverse momenta at RHIC and, in 
particular, LHC energies < px >= 0.39 — 0.61 GeV (see 
table [l] and [5]) the Jacobian has only a very small effect 
on the central source at sufficiently high values of \fs, and 
almost no effect on the fragmentation sources. 

For heavy-ion systems, the dependencies of the 
diffusion-model parameters on incident energy, mass and 
centrality at RHIC and LHC energies have been investi- 
gated in 19)^4 -26 . This Letter presents the first inves- 
tigation within the RDM for pp at high relativistic ener- 
gies. The corresponding parameters are shown in table [T] 
as functions of the cm. energy y/s. 

The time parameter Ti nt /T y is displayed as function of 
center-of-mass energy in table[TJ It is seen to increase with 
\fs. An increasing time parameter implies that the local 
maxima of the distribution function move further away 
from the beam rapidity (that increases with hi(y/s/m)) 
with increasing energy. In accordance with the expecta- 
tion, it indicates that the rapidity equilibration time t v 
decreases with rising energy, whereas the interaction time 
Ti n t depends only weakly on energy in the pp system. 




Fig. 1: Pseudorapidity distributions of produced charged 
hadrons in inelastic pp collisions at RHIC energies of ^/s — 
0.2 and 0.41 TeV as calculated with the three sources and fit- 
ted to PHOBOS data [6]. 



From the available data, it appears that the local max- 
ima occur at rather similar positions in pseudorapidity 
space, 77 ~ 2. This is characteristically different from 
heavy-ion collisions, where the maxima move outwards 
with increasing energy, as observed in AuAu data at RHIC 
[6j , and described in the RDM [9] . This qualitative differ- 
ence is most likely due to the larger spatial extent of the 
heavy-ion system. 

The partial widths (FWHM) as functions of energy are 
found to increase linearly with log y/s, table [l] Here the 
widths are effective values: beyond the statistical widths 
that can be calculated from a dissipation-fluctuation the- 
orem 27 within the RDM, they include the effect of col- 
lective expansion of the produced particles. The values 
at RHIC energies are resulting from a minimization with 
respect to the data that corresponds to the time evolution 
up to Tint- The integration is stopped at the optimum 
values of Ti n t/r v , Ti^.gg, and iV?? and hence, the explicit 
value of Ti n t is not needed. 

The normalization is given by the total number of pro- 
duced charged hadrons that is taken from experiment if 
available, or extrapolated in case of predictions at higher 
energies. Hence, the model contains four parameters. It 
provides an analytical framework to calculate the distri- 
bution functions, and to draw physical conclusions. 

The charged-particle distributions in rapidity space are 
obtained as incoherent superpositions of nonequilibrium 
and central ( "equilibrium" ) solutions of ([lj 



dN ch (y,t = Tint) 



dy 



= Nl h Rx{y,T in t) 



-N 2 ch R 2 {y, T int ) + N B c *R 3 (y, T mt ). (7) 



Results and discussion. — The results for pseudora- 
pidity distributions of produced charged hadrons in inelas- 
tic pp collisions at two RHIC energies of 0.2 and 0.41 TeV 
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Table 1: Three-sources parameters for pp collisions at RHIC energies (upper two lines) and at LHC energies (lower four lines). At 
RHIC energies the fragmentation sources from quark-gluon interactions with particle content dominate. At LHC energies 
the source from gluon-gluon collisions with particle content is the major origin of particle production at midrapidity. 
Midrapidity values (last column) are from PHOBOS (inelastic) |6 for |tj| < 1 at RHIC energies, and from CMS (NSD) [Ifl5] for 



77 1 < 0.5 at LHC energies. The 14 TeV value is calculated with the extrapolated parameters. See [H] for approximate average 
< Pt > —values. 
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Fig. 2: Pseudorapidity distributions of produced charged 
hadrons in inelastic pp collisions at RHIC energies of %/a =0.41 
TeV as calculated with the three sources and fitted to PHO- 
BOS data (6l. The underlying partial distribution functions 
including the Jacobian are shown as dash-dotted curves for 
the fragmentation sources arising from quark-quark and quark- 
gluon interactions, and as dashed curve for the central source 
arising from gluon-gluon interactions. The dotted curve does 
not include the Jacobian. 



Fig. 3: Charged-particle pseudorapidity densities in the central 
pseudorapidity region \ rj\ < 0.5 for non-single-diffractive (NSD) 
proton-proton collisions as function of the centre-of-mass en- 
ergy yfs. The value at 14 TeV (circle) is extrapolated from 
CMS data (dots, [5]). ALICE NSD data at 0.9 and 2.36 TeV 
are shown for comparison (triangles, Vm. Squares at RHIC 
energies of 0.2 and 0.41 TeV are inelastic PHOBOS pp— data 
for \q\ < 1 [6]. 



are shown in fig. [I] in comparison with PHOBOS data [6]. 
The three-sources model yields excellent agreement with 
the data. Here the overall normalization is taken from 
the data, and the fit parameters are the time parameter 
(that determines the mean values < y\$ >), the widths 
r ss , and the number of produced particles in the cen- 
tral source N?F. 

At RHIC energies, the multiplicity density at midrapid- 
ity has still a substantial contribution from the overlapping 



fragmentation sources. At 0.2 TeV, the contribution from 
the gluon-gluon source at 77 = is about 20 %, at 0.41 
TeV the midrapidity source is already much larger (48%), 
but the fragmentation sources still contribute 26% each, 
as shown in fig. [2] 

It should be mentioned that there exist detailed micro- 
scopic calculations of fragmentation sources from gq — > q 
and qg 



q diagrams by Szczurek et al. 



28j|29] for pion 
production in proton-proton and heavy-ion collisions at 
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Fig. 4: Pseudorapidity distributions of produced charged 
hadrons in pp collisions (NSD) at LHC cm. energies of 0.9, 
2.36, 7 and 14 TeV (bottom to top) as calculated in the three- 
sources approach and fitted to CMS NSD data [4][5]. At 0.9 
TeV UA5 NSD data are also shown [8], triangles. See fig. [5] for 
the underlying partial distribution functions at 7 TeV. 
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n 

Fig. 5: Pseudorapidity distributions of produced charged 
hadrons in pp— collisions at LHC cm. energy of 7 TeV. The un- 
derlying partial distribution functions including the Jacobian 
are shown as dash-dotted curves for the fragmentation sources 
arising from quark-quark and quark-gluon interactions, and as 
dashed curve for the central source arising from gluon-gluon 
interactions. The dotted curve does not include the Jacobian. 



SPS and RHIC energies. These processes are also respon- 
sible for the observed differences 30 in the production of 



positively and negatively charged hadrons, in particular, 
pions. An extension of these calculations to LHC energies 
is very desirable. 

To determine the RDM parameters in pp collisions at 
LHC energies, I have performed fits of the time parameters 
to the maxima of the double-humped charged-hadron dis- 
tributions, extrapolations of the partial widths T± t 2 gg with 
log \fs, and corresponding extrapolations of the number of 
produced particles in fragmentation and central sources as 
functions of log ^/s, see table[l] The number of particles in 
the central sources is at sufficiently high energy essentially 
determined by the measured pseudorapidity density near 
midrapidity that is plotted in fig. [3] as function of energy. 

The results at LHC energies are shown in fig. [4] The 
model results are compared with CMS data at 0.9, 2.36 
and 7 TeV [§[5], and UA5 data at 0.9 TeV (8). The calcu- 
lation at 14 TeV is performed based on an extrapolation 
of the multiplicity density at midrapidity with log \fs that 
yields dN/drj ~ 6.73 ± 0.30 at midrapidity. 

At LHC energies, the overall scenario changes in favor 
of particle production from the midrapidity source. The 
bulk of the midrapidity density is generated in the central 
source (73%) at 7 TeV, there is only a small overlap of the 
fragmentation sources at midrapidity as shown in fig. [4] 

In a comparison with calculations at LHC energies that 
do not include the Jacobian transformation as displayed by 
the dotted curve in fig. [5j it is evident that the midrapid- 
ity dip structure is essentially determined within the RDM 
by the interplay of the three sources for particle produc- 
tion, and only marginally influenced by the transformation 
from y— to 77— space at these high energies. The central 
distribution including the Jacobian has no dip at LHC en- 



ergies, but only a slight reduction in absolute magnitude 
at midrapidity, as shown by the dashed curve in fig. [5] 

There is, however, also the possibility that coherent 
QCD-effects contribute to the dip structure. Such effects 
go beyond the present calculation. They are also not con- 
sidered in numerical event generators (see [31] as an ex- 
ample), which provide rather accurate representations of 
RHIC data, in particular for transverse momentum dis- 
tributions. As compared to the analytical RDM these re- 
quire, however, a substantial numerical effort. 

Another purely empirical formulation of pseudorapidity 
distributions in multiple particle production at ^/s — 22.4 
to 1800 GeV based on several emitting centers along the 
rapidity axis had been given in [32] . It yields an analytical 
expression for the distribution function, and four param- 
eters are fitted to the data. As compared to the straight- 
forward physical interpretation of the three sources in 
the RDM it seems, however, difficult to assign a physi- 
cal meaning to the sources. 

The determination of the parameters within the RDM 
clearly goes beyond triple-gaussian fits that are modi- 
fied by the Jacobian, because the comparison with the 
data is based on, and constrained by, the underlying 
nonequilibrium-statistical description. Hence the depen- 
dence of the resulting parameters on incident energy as 
shown in table 1 is not arbitrary, but yields a consistent 
physical result. 

In particular, the time parameter increases with cen- 
ter of mass energy since the rapidity relaxation time de- 
creases. The width and particle content of the fragmen- 
tation sources do not change much with rising energy be- 
cause the number of contributing valence quarks stays con- 
stant, whereas the width and, in particular, the particle 



p-5 



G. Wolschin 



content of the central source that arises from gluon-gluon 
collisions increases substantially due to the large increase 
of gluons in the system at high energy and small values of 
Bjorken-x. 

As compared to the application of the RDM to heavy- 
ion collisions, it appears that transport phenomena are not 
fully developed in pp due to the small transverse size of 
the system. In the energy range from 0.2 to 14 TeV con- 
sidered here, the peak positions stay almost constant at 
r) ~ 2 in pseudorapidity space. The shift with energy that 
is present in AA systems 25 , and interpreted there as a 



multiparticle effect, does not seem to occur in pp. Hence, 
the full development of transport phenomena in higly rel- 
ativistic collisions requires a sufficiently large system in 
transverse size. 

Conclusion. — Based on the description of charged- 
hadron pseudorapidity distributions in pp collisions at 
RHIC and LHC energies in a nonequilibrium-statistical 
model, I have presented calculations of pseudorapidity dis- 
tributions of produced charged hadrons for pp collisions at 
RHIC energies of 0.2 and 0.41 TeV, and at LHC energies of 
0.9, 2.36, 7 and 14 TeV. These rely on the extrapolation of 
the transport parameters in the relativistic diffusion model 
(RDM) with increasing center-of-mass energy, and fits to 
the available data. 

In a three-sources model, the midrapidity source that is 
associated with gluon-gluon collisions accounts for about 
73% of the charged-particle multiplicity density measured 
by CMS at midrapidity in pp collisions at 7 TeV. The frag- 
mentation sources that correspond to particles that are 
mainly generated from valence quark - gluon interactions 
are centered at relatively large values of pseudorapidity 
(< Vi,2 >—< yi.2 >— T2.8) and hence, these contribute 
only marginally to the midrapidity yield. 

Since the Jacobian transformation from rapidity to 
pseudorapidity space is close to 1 at LHC energies 
due to the large mean transverse momenta, the size of 
the midrapidity-dip in the pseudorapidity distribution 
function is essentially determined by the relative particle 
content in the three sources, not by the Jacobian. Small 
corrections of the extrapolated values for the number of 
produced particles in the fragmentation sources may be 
required should measured distributions beyond pseudora- 
pidity rj — 2.5 become available from CMS, ATLAS and 
ALICE at LHC energies. 
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